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ABSTRACT 

We systematically examine how the presence in a binary affects the final core structure of a massive 
star and its consequences for the subsequent supernova explosion. Interactions with a companion star 
may change the final rate of rotation, the size of the helium core, the strength of carbon burning 
and the final iron core mass. Stars with initial masses larger than ~ HMq that experience core 
collapse will generally have smaller iron cores at the point of explosion if they lost their envelopes 
due to a binary interaction during or soon after core hydrogen burning. Stars below ~ HMg, on 
the other hand, can end up with larger helium and metal cores if they have a close companion, since 
the second dredge-up phase which reduces the helium core mass dramatically in single stars does not 
occur once the hydrogen envelope is lost. We find that the initially more massive stars in binary 
systems with masses in the range 8 — HM a are likely to undergo an electron-capture supernova, while 
single stars in the same mass range would end as ONeMg white dwarfs. We suggest that the core 
collapse in an electron-capture supernova (and possibly in the case of relatively small iron cores) leads 
to a prompt or fast explosion rather than a very slow, delayed neutrino-driven explosion and that this 
naturally produces neutron stars with low- velocity kicks. This leads to a dichotomous distribution of 
neutron star kicks, as inferred previously, where neutron stars in relatively close binaries attain low kick 
velocities. We illustrate the consequences of such a dichotomous kick scenario using binary population 
synthesis simulations and discuss its implications. This scenario has also important consequences for 
the minimum initial mass of a massive star that becomes a neutron star. For single stars the critical 
mass may be as high as 1O-12M0, while for close binaries, it may be as low as 6-8M Q . These 
critical masses depend on the treatment of convection, the amount of convective overshooting and the 
metallicity of the star and will generally be lower for larger amounts of convective overshooting and 
lower metallicity. 

Subject headings: binaries: close — globular clusters: general — supernovae: general — stars: neutron 
— X-rays: stars 



1. INTRODUCTION 

For the last decade, it has generally been accepted 
that the high speeds inferred for isolated, normal ra- 
dio pulsars imply that neutron stars (NSs) receive a 
large impulse, or "kick," at birth. Measured proper 
motions for ~100 radio pulsars indic ate typical kick 
speed s in excess of 100 -200 km s -1 (|Lyjie_^JLiQrim£rl 
lltlMlmmsen fc Phinnevlll997tlCordes fc Chernofi1ll99St 
lArzoumanian. Chernoff. &; Cordesl l2002jh though the 
functional form of the underlying speed distribution is 
poorly constrained. The physical mechanism that causes 
this kick is presently unknown, but is presumably the 
result of some asymmetry in the core colla pse or sub- 
sequen t supernova (SN) explosion (see, e.g.. lPfahl et alJ 
( 20023) for discussion and references). 

In apparent conflict with the high speeds of isolated 
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radio pulsars, iPfahl et alJ l)2002c|) identified a new class 
of high-mass X-ray binaries (HMXBs) wherein the neu- 
tron stars must have been born with relatively low kick 
speeds. These HMXBs are distinguished by their long 
orbital periods (P or b > 30 d) and low eccentricities (e < 
0.2), with the prime example being X Per/4U 0352+309 
(P or b — 250 d; e ~ 0.1). The orbits of these systems 
are sufficiently wide that tidal circularization is negligi- 
ble, so that the observed eccentricities should reflect the 
conditions immediately after the SN explosion. Such low 
eccentricities essentially require that the neutron stars in 
these systems received low kick speeds of <50kms _1 . 

Further observational evidence that at least some neu- 
tron stars receive low kicks at birth is provided by the 
fact that a large number of neutron stars are found 
in globular clusters, where some massive globular clus- 
ters may contain more than ~ 1000 ne utron stars 
(Pfahl. R appaport. fc Podsiadlowskil l2002b|L Since the 
central escape velocity is generally < 50kms _1 , essen- 
tially all of the neutron stars born in a globular cluster 
should escape from the cluster if they received a kick 
consistent with the kick distribution for single radio pul- 
sars (for a detailed discussion of this so-called 'neutron- 
star retention problem' see Pfahl et al. 2002b). If all of 
these neutron stars were originally born in massive bi- 
naries, this would alleviate the problem somewhat, since 
in this case the momentum imparted to the neutron star 
would be shared by the total mass of the system, lead- 
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ing to a corresponding ly smaller space velocity of th e 
post-supernova binary ijBrandt fc Podsiadlowskil il995'). 
However, while the membership in a binary increases the 
number of neutron stars that can remain bound in the 
cluster significantly, the effect may not be large enough to 
explain the observed numbers unless clusters were much 
more massive at an earlier epoch than they are today 
(Drukier 1996; Pfahl et al. 2002b). 

To simultaneously account for the new class of H MXBs 
and the high speeds of radio pulsars. [Pfahl et al] {2002c) 
suggested that neutron stars originating from progeni- 
tors that are single or members of wide binaries receive 
the conventional large kicks, while neutron stars born 
in close binaries receive small kicks (also see Katz 1975, 
1983 and Hartman 1997 for earlier more ad hoc sugges- 
tions of a significan t population of neu tron stars with 
small natal kicks). iPfahl et all l)2002c|) further argued 
that the proposed dichotomy between high and low kick 
speeds, and its relation to the evolutionary history of the 
NS progenitor, might be associated with the rotation rate 
of the collapsing core. However, as we will show in this 
paper, the core structure itself (e.g., its mass and com- 
position) depends strongly on whether a star evolves in a 
close binary or in isolation. This can produce differences 
in the actual supernova and may allow a prompt (or at 
least a fast) supernova explosion mechanism to be suc- 
cessful in some cases; this in turn may lead to relatively 
low supernova kicks. 

The outline of this paper is the following: In § 2 we 
systematically discuss the differences in the core evolu- 
tion for stars in close binaries and in single systems/ wide 
binaries and discuss the implications for the core-collapse 
phase. In § 3 we develop a general scheme for forming 
neutron stars in different types of systems and the ex- 
pected differences in kick velocity, associating them with 
individual known systems or classes of systems. In § 4 
we discuss how this scheme can be tested both observa- 
tionally and suggest how further progress can be made 
on the theoretical side. 

2. BINARY EVOLUTION AND THE PRE-CORE-COLLAPSE 
CORE STRUCTURE OF MASSIVE STARS 

It is often naively assumed that the evolution of he- 
lium cores is the same irrespective of whether the core is 
surrounded by a hydrogen envelope or not, and that the 
final core structure will be similar in the two cases. How- 
ever, binary evolution may affect the final pre-supernova 
structure of a massive star in several fundamentally dif- 
ferent ways: (1) the rate of rotation of the immedi- 
ate pre-supernova core, (2) the size of the helium core, 

(3) the occurrence of a second dredge-up phase at the 
beginning of the asymptotic-giant branch (AGB), and 

(4) the C/O ratio at the end of helium burning which 
affects the strength of carbon burning and the final 
size of the iron core. These effects can dramatically 
change the condition of core collap se as first pointed 
out bv i B rowm L ee, fc Bethel (^999) and in particular 
iBrown et alJ 1)20011) . who showed that the final iron cores 
of massive stars will be significantly smaller for stars that 
have lost their hydrogen-rich envelopes soon after the end 
of core hydrogen burning. 

The role of rotation 
While massive stars are generally rapid rotators on the 



main sequence, there are several efficient mechanisms by 
which they lose their angular momentum during their 
evolution. The final rotation rate of the core of a mas- 
sive star depends on whether the star passed through a 
red-supergiant phase during which the core will be effec- 
tively braked by the hydrodynamic and magnetic cou- 
pling to the slowly rota ting envelope ijSoruit fc Phinnevl 
119981: iHeeer et alJl2003^ and the wind mass loss during 
a subsequent Wolf-Rayet/helium-star phase, which can 
be very effic ient in extracting angula r momentum from a 
helium star llHeeer fc Wooslcv 2003) and slowing it down 
in the process. IPfahl et al.l l)2002oD argued that stars that 
lost their H-rich envelopes soon after the main-sequence 
phase (so-called early Case B mass transfer 6 ) might avoid 
this phase where the core is effectively braked and may 
still be rapidly rotating at the time of core collapse. 

However, whether the star can avoid spin-down also 
depends on whether mass transfer is dynamically stable 
or u nstable. In the case of stable Ca s e B Ro che-lobe over- 
flow, Eiflii^^^nsiiI^^^£iI^xI3 <|2003al) showed that, 
if the mass-losing star remains tidally locked to the or- 
bit, this provides a very efficient method of slowing down 
the rotation rate of the mass-losing star. On the other 
hand, in the case of unstable mass tra nsfer, leading to a 
common-envelope and spiral- in phase ijPaczvriskll 1 976j) . 
the spiral-in timescale is much shorter than any realistic 
synchronization timescale, and tidal locking would not be 
expected. In this case, the core of the mass-losing star 
could still be rapidly rotating after the ejection of the 
common envelope. This imp lies that the scenario sug- 
gested bv lPfahl et alJ (|2002c) probably requires late case 
B mass transfer (i.e., dynamically unstable mass transfer 
that leads to a CE phase before the core has been spun 
down significantly). 

Another situation which may lead to a rapidly ro- 
tating core is the complete slow merger of two mas- 
sive stars, in particular if it occurs after helium core 
burning, as in models for the p rogenitor of SN 1987A 
ijlvanova fc Podsiadlowskil |2003) which predict a very 
rapidly rotating core for the immediate pre-supernova 
star (also see Joss & Becker 2003). 

Finally we note that, if the exploding star is still accret- 
ing from a companion at the time of the s upernova, one 
would also expected a rapidly rotating core i|Langer et alJ 
2003). 

As this discussion shows, binary interactions can sig- 
nificantly affect the final pre-supernova rotation rate of 
the core of a massive star, although the details can be 
rather involved and are not completely understood. In 
addition, what is even less clear at the present time is 
how this affects the physics of the core collapse itself 
and, in particular, the magnitude of the kick imparted 
to the newborn neutron star. 

The size of the helium core and the second dredge-up 

The final size of the helium core depends strongly 
on the evolution of the H-rich envelope. During core 

6 It is common practice to distinguish among three evolution- 
ary p hases of the primary at the o nset of mass transfer, fol- 
lowing KiDDcnhahn & Wciecrt (1967) (see also Lauterborn 1970; 
Podsiadlowski. Joss. & Hsu 1992:). Case A evolution corresponds 
to core hydrogen burning, Case B refers to the shell hydrogen- 
burning phase, but prior to central helium ignition, and Case C 
evolution begins after helium has been depleted in the core. 
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helium burning, the helium core ordinarily grows sub- 
stantially because of hydrogen burning in a shell - of- 
ten the dominant nuclear burning source - which adds 
fresh helium to the core. On the other hand, in a 
binary a star may lose its H-rich envelope before He 
burning (or early during the He-burning phase). In 
this case, the helium core can no longer grow and 
may in fact shrink because of the strong stellar wind 
expected in the subsequent Wol f -Rayet phase (e.g., 
Wooslev . Langer. fc Weaver! Il995a t IWellstein fc Langer 
19991 iWellstein. Langer. fc Braunl I200H iPols fc Dewi 



2002). Therefore, the final mass of the helium core will 



often be lower for stars in close binaries than in single 
systems/wide binaries (see Fig. 1). 

Another factor that strongly affects the mass of the fi- 
nal helium core is the occurrence of a second dredge-up 
phase. When stars up to ~ 11M Q ascend the asymp- 
totic giant branch, they generally experience a second 
dredge-up phase where the convective envelope pene- 
trates deep into the H-exhausted core and dredge s up 
a significant fraction of the helium core l)Ibenlll974j) . As 
a consequence, the size of the helium core can be dra- 
matically reduced (by up to ~ 1.6 M Q ; see Fig. 1). How- 
ever, the occurrence of the second dredge-up depends 
on the presence of a convective H-rich envelope. If the 
star loses its H-rich envelope before this phase, dredge- 
up does not occur. In this case, the final size of the 
helium core is larger for a star in a binary than its sin- 
gle counterpart (note that this is the opposite of what 
happens to their more massive counterparts). This is il- 
lustrated in Figure 1 which shows the final helium core 
mass as a function of initial main-sequence mass for sin- 
gle stars (thick solid curve; iPoelarends fc Langeill200J) 
and for stars in binaries that lose their envelopes either in 
Case A or Case B mass transfer ( based on the results of 
IWellstein. Langer. fc Braunl 1200 ll) . The almost discon- 
tinuous change of the final helium mass around 12 Mq 
is a direct consequence of the fact that stars below this 
mass have experienced a second dredge-up phase, while 
stars above this mass do not or only dredge up a moder- 
ate amount of the helium core. Note that the final mass 
for these stars is less than the minimum mass for core col- 
lapse (~ 1.4 M Q ). After the dredge-up phase, the helium 
core may grow again because of hydrogen shell burning, 
just as in a normal AGB star. Whether the core can 
reach the critical mass for core collapse depends on the 
timcscale on which the star loses its envelope in a stellar 
wind or superwind. While this is somewhat uncertain, 
we estimate that for the models in Figure 1, single stars 
below ~ 12 Mq will not reach the critical mass for core 
collapse and will end their evolution as ONeMg white 
dwarfs. There may be a small mass range around this 
boundary where single stars reach the condition for core 
collapse. 

In contrast, the shaded region between the dot-dashed 
thick lines indicates the expected range of the final he- 
lium core mass for stars that lose their envelopes by bi- 
nary interactions (in this case, the final core mass de- 
pends on the evolutionary phase and the core mass at 
the time of m ass trans f er). 

In the past. lNomotol ((19841 IT987kJ| has argued that an 
electron-capture supernova is the expected fate for stars 
with main-sequence masses in the range of 8-IOMq. In 
his calculations, these stars developed helium cores in the 
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Fig. 1. — Final mass (thick solid line) and maximum mass 
(thick dashed line) of the he lium core in single sta r s as a function 
of initial mass according to Poclarcnds & Langer (2004) (PL04), 
extrapolated for initial masses above 12.5 Mq. The final helium 
core masses from the calculations of Wooslev & Weaver ( 1995) 
are indicated by a thin solid line. The final helium core masses 
of close binary models undergoin g Case B mass transfer from 
IWellstein. Langer. fc B raun (2001) (WLB01) are shown as the up- 
per dot-dashed line, while those experiencing Case A mass transfer 
(WLB01) may lie anywhere between the lower dot-dashed line and 
the Case B line. The results from the binary calculations have been 
extrapolated for initial masses below 12 Mq (dotted part). Note 
that the PL04 and WLB01 models have been computed with the 
same assumptions for convective mixing, while the WW95 mod- 
els assumed a higher semiconvective mixing efficiency. The light 
dashed horizontal lines give the range for the final helium core mass 
for which the star may experience an electron-capture supernova. 
Note that the parameter range for which this may occur for a single 
star is very small. 



range of M — 2.0-2.5 Mq 7 and never developed an iron 
core; in this case, the collap se is triggered by the e lectron 
capture on 24 Mg and 20 Ne ilNomotdll98l Il987al) (for a 
recent discussion see lWanaio et al.M2003|) . In Figure 1, 
the light dashed horizontal lines indicate approximately 
the range of helium-core masses which can be expected 
to lead t o an electron-capture s upernova (based on th e 
results of lNomotd (119841198731) : also see lHabetslll986|h 
As Figure 1 shows, because of the second dredge-up, the 

7 We note that in Nomoto's calculations, stars in the range 
of 8-IOMq either did not experience a second dredge- up phase 
or only during carbon shell burning. As a consequence, most 
of his models in this mass range, unlike ours, experienced a 
core-collapse supernova even in the case of single stars. We 
suspect that this difference can be attributed to the different 
opacities employed. The new OPAL opacities are significantly 
larger in the critical temperature range, which makes dredge-up 
more likely or occur earlier. We emphasize that this dredge-up be- 
haviour is qualitatively found in al l other recent, detailed studies 
of s tars in thi s mass range iRi tossa. Garcfa-Berro. fc Ibenl 
[1993; ICiarcia-i icrro. Kitossa. fc lb en| U997J; 

llben. Ritossa. fc Garcfa-Berrol 119971 lEldridge fc TouH 120041 ). 



We note, however, that the range of initial masses which leads to 
helium cores in the range of 2.0- 2.5 Mq depends on some of the 
uncertainties in the stellar modelling, in particular the treatment 
of convection and convective overshooting (for further discussion 
of the uncertainties see § 2.2). 
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mass range for which single stars experience an electron- 
capture supernova may be very narrow if non-existent, 
while there is a large mass range for which it may occur 
for a star in a close binary. Indeed a binary channel may 
be the only place where it can occur. 

The C/ O ratio and the strength of carbon burning 

Another more subtle effect is that the lack of a H- 
burning shell leads to a lower C/O ratio at the end of 
helium core burning which affects the strength of sub- 
sequent carbon burning and the final size of the iron 
core. This ef f ect w as first pointed out and explained by 
iBrown et alJ l)2001l) and can be understood as follows. 

In the late phase of helium core burning, the alpha- 
capture reaction C 12 +a — > O 16 , which destroys carbon, 
tends to dominate over the carbon-producing triple-a re- 
action (because of the different functional dependence on 
the number density of a-particles). As a consequence, 
carbon is systematically being destroyed in this late core- 
hclium-burning phase (this switch occurs when the he- 
lium mass fraction falls below Y = 0.1-0.2). In he- 
lium cores surrounded by a hydrogen-burning shell, both 
the hydrogen-exhausted core and the convective helium- 
burning core continue to increase during helium core 
burning. This leads to the continued injection of fresh 
helium into the helium-burning core and prolongs the 
phase in which carbon is preferentially destroyed. This 
eventually produces a smaller carbon abundance at the 
end of core helium burning than would be the case for 
a naked helium core, where the convective core does not 
grow. This has drastic consequences for the subsequent 
carbon-burning phase if helium or carbon-oxygen cores 
of similar sizes are compared. For high carbon abun- 
dances (as expected for naked helium cores), the phase 
of convective carbon shell burning lasts longer, typically 
leading to smaller carbon-exhausted cores). This in turn 
produces smaller iron cor es with steeper density gradi- 
ents o utside the iron core l)Frver et alJl2002t iHeeer et al.1 
120021) . Thus, for the same initial size of the helium or 
carbon core, the higher carbon abundance in a star that 
lost its hydrogen-rich envelope before the central helium 
abundance dropped below ~ 10 % will result in a pre- 
supcrnova structure that more easily produces a success- 
f ul supernova. 

IBrown et alJ l|2001j) demonstrated that this dichotomy 
leads to much smaller iron cores for massive stars and 
that, as a consequence, even a 60 Mq main-sequence star 
may produce a neutron star rather than a black hole if 
it has lost its envelope before its helium core-burning 
phase. In contrast, the minimum mass for blac k-hole 
formation for single stars may be as low as 20 Mq l|Frverl 
1999; iFrver fc Kalogerall200l(l . 

In analogy to these results, it is reasonable to expect 
that there will be significant differences in the core prop- 
erties even for massive stars with M < 20 Mq, possibly 
allowing for s uccessful prompt (fast ) supernova explo- 
sions (see e.g., (Sumivoshi et al. 200l]and § 2.1). 

2.1. Prompt (fast) explosions and supernova kicks 

At present, neither the mechanism that produces a suc- 
cessful core-collapse supernova nor t he origin of super- 
nova kicks is properly understood (see lJanka et alJl2003t 
lFrverfl2003l for detailed recent reviews, and also Fryer & 



Warren (2002, 2003)). In one of the most popular ex- 
plosion scenarios, it is the delayed heating by neutrinos 
that revives the outgoing supernova shock several 100 ms 
after it stalled in the initial core bounce (i.e., after many 
dynamical timescales). In this scenario, the origin of the 
supernova kick may be connected with asymmetries in 
this long phase where the explosion develops (e.g., due 
to the continued accretion onto the proto-neutron star or 
caused by t he strong convect ion in the neutrino-heated 
region; e.g., Uankaet al.1 12003'). In contrast, in a prompt 
supernova explosion, the initial bounce drives a success- 
ful supernova shock on the dynamical timescale of the 
proto-neutron star. The absence of a long phase where 
the explosion teeters at the brink of success could then 
be the cause for the absence of a large supernova kick. 

Such a scenario is supported by some of the most recent 
core-collapse simulations by Scheck et al. (2004) (see, 
however, also Fryer & Warren 2003 for a very different 
view) . In the simulations by Scheck et al. (2004) , a large 
supernova kick (up to and above 500 km s _ ) is caused by 
asymmetries in the neutrino flux, which have their origin 
in low-order instabilities driven by the convective motion 
behind the stalled shock. An essential requirement in 
these simulations, which allows these convective insta- 
bilities to grow, is that the duration of the convective, 
stalled phase is longer than ~ 500 ms (i.e., many convec- 
tive turnover timescales). Such slow explosions are ex- 
pected for fairly large iron cores. In contrast, for a small 
iron cores or in the case of electron-capture supernovae, 
the simulations by Scheck et al. (2004) suggest fast ex- 
plosions where these convectively driven instabilities are 
unable to grow, leading to rather small kick velocities 8 . 

It seems quite attractive to relate the dichotomy in the 
supernova kicks to the differences between a slow and a 
prompt (fast) supernova explosion. Whether a supernova 
explosion develops promptly or in a delayed manner de- 
pends mainly on the difference between the mass of the 
Fe-Ni core and the mass of the collapsing core (which 
in turn also depends on the initial entropy in the core), 
since this determines the amount of shock energy that is 
cons umed in the nuclear dissociati o n of h eavy elements 
(see Iffillebrandt. Nomoto . &: Wolfl ill 9841) for a detailed 
discussion and references; iSumivoshi et al.ll2001|) . 

Electron-capture supernovae provide a particularly 
promising scenario for a prompt (fast) explosion, since 
the whole core collapses to nuclear densities; this makes 
it much easier for the shock to eject the envelope, pre- 
venting the growth of the instabilities that lead to large 
kicks. As our previous discussion shows, electron-capture 
supernovae may only (or mainly) occur in close binaries; 
in this case, neutron stars with low kicks may be (al- 
most) exclusively produced in close binary systems (with 
orbital periods < a few 100 d). 

2.2. The minimum mass for core collapse 

An important related issue is the question of the min- 
imum initial mass, M m j n , of a star that leads to a core- 
collapse supernova in a single or binary system. It is 
commonly assumed that this minimum mass is around 

8 Note that, unlike the case of a prompt explosion, the fast explo- 
sions in the simulations by Scheck et al. (2004) occur on a timescale 
long compared to the dynamical timescale of the proto-neutron 
star. 
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8 Mq, the minimum mass above which ignite carbon off- 
center (rather th an explosiv ely in the center) and form 
an ONeMg core i|Ibenlll974|) . If the ONeMg core is able 
to grow to reach the Chandrasekhar mass, it will col- 
lapse in an electron-capture supernova. However, if a 
star experiences mass-transfer already on the main se- 
quence, an initial star as massive as ~ 16 Mq may end 
its evolution as a white dwarf rather than experience 
core collapse 9 . On the other hand, as discussed above, 
stars below 12 Mq may only experience core collapse if 
they have lost their envelopes by binary interactions af- 
ter their main-sequence phase but before experiencing 
dredge-up on the AGB. This implies that the value of 
M m - ln may be as high as 12 Mq for single stars and as 
low as 8 Mq for relatively close binaries. 

The exact value for M m j n is quite sensitive to the 
treatment of convection and, in particular, the amount 
of convective overshooting, and the metallicity of the 
star. The value of M m - m = 12 Mq was obtained for 
models that used the Ledoux criterion for convection 
without convective overshooting. We estimate that us- 
ing the Schwarzschild criterion would reduce M m i„ to 
~ 11 -Mq if no convective overshooting is included, and 
to ~ 10 Mq if the recent empirical calibration of convec- 
tive overshooting (jPols^t^aJJIWgjjjS^ 1997) 
is adopted (also see | Ritossa. Garcia-Berro. fc Iber]ll999t 
Eldr idge fc foul I2004J) . Similarly, the minimum mass 
for off-center carbon ignition may be as low as ~ 6 Mq 
for models including convective overshooting (Han 2002, 
u npublished). 

lHan. Podsiadlowski. fc Eggletonl l)1994|) also found 
that, for low-metallicity [Z = 0.001) models without con- 
vective overshooting, the minimum mass for off-center 
carbon ignition was systematically lower (~ 6 Mq) than 
for solar metallicity (~ 8Mq). 

As these discussions illustrate, the initial mass range 
that leads to the formation of ONeMg white dwarfs, 
which also determines the minimum mass for stars that 
will experience core collapse, depends both on the details 
of binary interactions and on the stellar properties. Dif- 
ferent treatments of convection are expected to lead to 
initial mass ranges ranging from [6,9] Mq to [9,12] Mq. 
In addition at low metallicity, these ranges should be 
systematically shifted t owards lower masses (by perhaps 
2M m for Z = 0.001; lHan. Podsiadlowski. fc Eggletonl 
1994). Considering that this is an important parame- 
ter in galactic modelling, a thorough re-examination of 
this issue is urgently n eeded ijPoelarends fc Langerll2004l 
lEldridge fc Tout1l200l . 

We emphasize that the exact values of this mass range 
do not affect the main arguments in this paper, since this 
only shifts the mass range in which an electron-capture 
supernova can be expected, but does not change the ex- 
pected dichotomous behaviour. 

3. A DICHOTOMOUS KICK SCHEME 

The scenario for NS kicks proposed herein has a signif- 
icant impact on the theoretical production probabilities 

9 This occurs when the system experiences an ad- 
ditional RLOF phase a fter core He exhaustion (Case 
ABB mass transfer: IX Vellstein, Lang er. fe Braunl 120011 : 
[Podsiadlowski, Rappaport, Si Han 2003) during a He red- 
supergiant (RSG) phase; strong mass loss during the He-RSG 
stage may lead to a similar outcome. 
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Fig. 2. — Results of a single HMXB binary population synthesis 
calculation that utilizes our proposed dichotomous kick scenario. 
Each red (blue) dot represents a system where the primary was 
(was not) highly evolved when it lost its hydrogen-rich envelope, 
and the exposed core evolved to collapse to form a neutron star 
with a subsequent large, conventional (small, unconventional) na- 
tal kick; see the text for details. Markers in the right panel indicate 
the observed wide, low-eccentricity HMXBs {filled circles), and the 
well-known eccentric HMXBs {filled triangles). The filled squares 
show the three radio pulsars with massive binary companions in ec- 
centric orbits. We did not include observed systems with -P rb <10 
days, since their orbital parameters arc likely to have been altered 
by tidal circularization effects. 



and distributions of orbital parameters of binary systems 
containing neutron stars. We illustrate this by means of 
a Monte Carlo binary population synthesis (BPS) cal- 
culation. Below we provide a brief description of the 
importan t elements of the code; further deta i ls may b e 
found in Pfahl . Rappaport. fc Podsiadlowskil (I2002albf): 
Pfahl et alJ l|2002cj) ; IPfahl. Rappaport. fc Podsiadlowskil 
(2003). The BPS code follows a randomly generated 
sample of massive primordial binaries through the phase 
of mass transfer from the primary to the secondary. 10 
The initial primary and secondary masses, Mu and M^i, 
are drawn from the respective distributions p{Mu) cx 
M^ 2 ' 5 and p(q l ) = 1, where q l = M 2 i/M u < 1 is the 
initial mass ratio. For simplicity, we assume circular or- 
bits, and the initial orbital separation, ctj, is chosen from 
p(ai) oc ar 1 . 

Given some critical mass ratio, which we take to be 
q c = 0.5, the mass transfer from the primary to the sec- 
ondary is assumed to be stable if > q c and the envelope 
of the primary is mostly radiative when mass transfer be- 
gins (so-called Cases "A" , "early B" , and "early C" ) , and 
dynamically unstable if qi < q c or the primary has a con- 
vective envelope (Cases "late B" and "late C"). In any 
case, we assume that the entire hydrogen-rich envelope 
of the primary is removed during mass transfer, leaving 

10 Here "primary" and "secondary" refer to the initially more 
and less massive star, respectively. 
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Fig. 3. — Illustrative binary evolution scenarios leading to the 
formation of high-mass X-ray binaries. In the left panel, the binary 
is sufficiently wide at the onset of mass transfer that the primary 
has a fully formed 3.9 Mq He core. By contrast, in the right panel 
the primary's envelope is lost while the He core is substantially 
lower in mass (2.4 Mq). In the proposed dichotomous kick scenario 
the former case leads to a large natal kick, while the latter case 
results in a "prompt (fast)" supernova event and a smaller kick 
(see text). 



only the primary's hydrogen-depleted core. 

For stable mass transfer, we suppose that the sec- 
ondary accretes a fraction, (3 = 0.75, of material donated 
by the primary, and that the remaining mass escapes the 
system with a specific angular momentum that is a = 1.5 
times the orbital angular momentum per unit reduced 
mass. The value of a = 1.5 is characteristic of mass 
loss through the L2 point, but the choice of f3 = 0.75 
is fairly arbitrary - although it is convenient, as it gives 
af/cii ~ 1, where a/ is the final orbital separation. Dy- 
namically unstable mass transfer results in a common- 
envelope phase, wherein the secondary spirals into the 
envelope of the primary. The common-envelope phase 
was treated in the same way as in Pfahl et al. (2002c). 
A fraction, tjce is 1, of the initial orbital energy is avail- 
able to unbind the common envelope from the system. If 
insufficient energy is available, the two stars will merge. 
Otherwise, the envelope of the primary is dispersed, the 
secondary emerges near the ZAMS without having ac- 
creted any significant amount of mass, and the orbital 
separation is ^100 times smaller. A merger results in 
nearly every case where < q c and the primary's enve- 
lope is radiative when mass transfer starts. A number 
of these evolutionary steps are illustrated in the top six 
panels of the schematic shown in Figure 

Upon exhausting its remaining nuclear fuel, the ex- 
posed core of the primary explodes as a Type lb or Ic SN. 
We assume that the supernova mass loss and natal NS 
kick are instantaneous, and that the orientations of the 
kicks are distributed isotropically. In light of the discus- 



sion above, the specific prescription adopted for NS kicks 
is as follows. Kick speeds are drawn from a Maxwellian 
distribution, with a mean that depends on the initial 
mass of the primary and its evolutionary state when it 
first fills its Roche lobe (see Fig. El- If 8 < M u < 14M 
and mass transfer begins before helium ignition in the 
primary's core (Cases A and B), the mean kick speed 
takes a small value of 30 km s _1 , while in all other cases 
the mean is 300 km s _1 . The range of masses and the 
actual mean kick speeds used here are reasonable choices, 
but are somewhat arbitrary, and were chosen mainly for 
illustrative purposes. 

The results of our BPS calculation are shown in Fig- 
ure [5] in the P or b — M2 and P or b — e planes, where M2 
is the mass of the secondary star, including the mass it 
has accreted from the primary, and e is the orbital ec- 
centricity. Each small dot represents an incipient high- 
mass X-ray binary immediately after the supernova ex- 
plosion that produces the neutron star. The red dots 
represent systems which experienced the nominal natal 
neutron star kicks (see also left column of Figure|3J) while 
the blue dots are systems where the primary lost its en- 
velope early enough so that the natal kicks were much 
smaller - in accord with the dichotomous kick scenario 
proposed herein (see right hand panel of Fig. EJ. 

The vertical blue strip in the P or b — M% plane of Fig. 
[21 results mostly from early Case B systems in which the 
primaries had largely radiative envelopes and start mass 
transfer before core He ignition. The resultant mass 
transfer is stable and the H-exhausted core evolves to 
a "prompt (fast)" SN explosion with a small natal kick. 
These systems generally acquire small eccentricities, es- 
pecially for P or b in the range of 10-100 days. These 
are presumably the systems that evolve to become the 
newly identified class of HMXBs with relatively wide, 
low-eccentricity orbits (e.g., X Per/4U 0352 + 30; 7 
Cas/MX 0053+604; XTE J0543-568; the filled circles in 
Fig. 2). For the systems of this type with still wider 
orbits, even the small natal kicks assumed here are suffi- 
cient to induce substantial eccentricities. 

The blue "tail" -shaped region in the P or b — M2 plane 
of Figure [5|results from late Case B systems commencing 
mass transfer when the primary had a large, convective 
envelope, but still before core He core ignition. These 
lead to common-envelope events and result in systems 
with short orbital periods (i.e., P or b < 10 days; see e.g., 
Fig. |2J) . Since such short period systems would, in any 
case, be circularized by tidal friction, it may be difficult 
to distinguish these from their Case C cousins where the 
primaries were much more evolved at the start of mass 
transfer, and which received larger natal kicks. Systems 
in these latter categories (late Case B and late Case C) 
presumably include many of the famous HMXBs such 
as Cen X-3, LMC X-4, 4U 0900-40, and SMC X-l, with 
current-epoch (approximately) circular orbits and short 
orbital periods. The red dots in Figurc[2]are those receiv- 
ing a large kick, where the primary was highly evolved 
at the onset of mass loss (Case C). These likely account 
for the class of eccentric HMXB systems with P or b in 
the range of - 20 - 200 days, e.g., 4U 0115+63, GX 
301-2 , 4U 1145-619, and EXO 2030+375 (|Bildsten et alJ 
Il997|) (filled triangles in Fig. 2) and the three radio pul- 
sars with massive companions in highly eccentric systems 
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(filled squares in Fig. 2) PSR 1259-63 (Johnston et al. 
1992), PSR J0045-7319 (Kaspi et al. 1994, 1996) and 
PSR J1740-3052 (Stairs et al. 2001). 

Finally, we note that this dichotomous kick sce- 
nario helps the long-standing "retention problem" for 
neutr on stars in globular clusters. In an earlier 
work ( Pfa hl. Rappaport. fc Podsiadlowski 2002^) we re- 
ported that our original dichotomous kick scenario 
yielded an enhancement factor of about 4 over the 
fraction of retained neutron stars without the appro- 
priate subset of small kicks. We also pointed out 
(Pfahl. Rappa port. fc Podsiadlowski l2002bh that a di- 
chotomous kick scenario would likely increase the theo- 
retical formation rate of double neutron star systems by 
approximately an order of magnitude. Both of these en- 
hancements are expected to be carried over to this newer 
kick scenario. 

4. DISCUSSION AND FUTURE WORK 

As we have shown, the presence in a binary can dra- 
matically affect the structure of the core of a massive star 
at the time of core collapse. Stars above ~ 11M are 
generally expected to have smaller iron cores if they lose 
their envelopes in a close binaries. Stars in the range of 
8-11 Mq may explode in an electron-capture supernova 
if they are in a close binary, while single stars or stars in 
wide binaries will experience a second dredge-up phase 
and are more likely to end their evolution as ONeMg 
white dwarfs. We suggested that in the case of small 
iron cores or in the case of an electron-capture super- 
nova, the supernova occurs through a prompt (fast) ex- 
plosion rather than a delayed neutrino-driven explosion 
and argued that this is more likely to produce neutron 
stars with low kick velocities. 

While speculative at the moment, this scenario has 
important observational implications and suggests the 
need for futher theoretical studies. These include a sys- 
tematic exploration of the late evolution of the cores 
of stars around 10 Mq, both for stars evolved in isola- 
tion and in a close binary, the dependence on metal- 
licitv and the amoun t of convective oversho oting (e.g., 
iPoelarends fc Langerll200l lEldridge fc ToutJl2Tj04|) . and 
a re-examination of the core-collapse phase for electron- 
capture supcrnovae and low-mass iron cores. One im- 
mediate prediction of this scenario is that neutron stars 
that formed from ONeMg or low-mass iron cores should 
produce neutron stars of systematically l ower mass (~ 
1.25 M m for a typical equation of state 11 ). iWanaio et alJ 
(2003) suggested that prompt explosions are the site 
for the r-process; if some of this r-processed matter is 
captured by the companion stars, this may produce de- 
tectable chemical anomalies in neutron-star companions 
(e.g., X Per) similar to the case of the companion in 

The pulsar that formed last in the recently discovered double- 
pulsar binary PSR J0737-3039 [Burgay et al. 2003; Lyne et al. 2004] 
has a low mass of ~ 1.25 Mq, although at present it is not clear 



Nova Scorpii ijlsraelian et al.lfl999t iPodsiadlowski. et alJ 
2002). Note also that, in this case, only small amounts 
of iron and other heavy elements are ejected in the su- 
pernova. If all neutron stars in globular clusters were 
to form through such a channel, one would not expect 
significant chemical pollution of the cluster with heavy 
eleme nts from supernova ejecta (jPrice fc Podsiadlowskil 
1993). 

Our scenario also has important implications for the 
minimum mass for neutron-star formation; it suggests 
that only single stars more massive than ~ 10-12 Mq 
become neutron stars, while in binaries the mass may be 
as low as 6 - 8 Mq . The exact values will depend both on 
the amount of convective overshooting and the mctallic- 
ity of the population; larger amounts of overshooting and 
lower metallicity are both expected to shift these critical 
masses to lower values. The detection of both young, 
single massive white dwarfs and neutron stars in bina- 
ries in open clusters with a turnoff mass around 10 Mq 
would provide direct evidence for such a dichotomy. An- 
other promising method to observationally constrain the 
minimum mass for a core-collapse supernova is through 
the detection of the progenitors of Type II- P supernovae, 
stars that have not lost all of their H-rich en velopes at the 
time of the supernova explosion (see e.g.. ISmartt et aD 

Finally, it is worth mentioning another channel which 
is likely to produce an electron-capture supernova and 
possibly a neutron star with a low kick if the ba- 
sic picture described in this paper is correct: this in- 
volves the accretion-induced collapse (AIC) of a white 
dwarf to form a neutron star, either because an ONeMg 
white dwarf accreting from a companio n star is pushed 
above the Chandraskekh ar limiting mass l)Nomotol l987b: 
INomoto fc KonddH99l or as the result of the merger of 
two CO white dwarfs which is also likely to lead to the 
formatio n of a neutron star ra ther than a Type la su- 
pernova l)Nomoto fc Ibeniri985j) . In principle, AIC could 
produce single low-vel ocity neutron stars with a rate a s 



high as 3 x lO^yr" 1 l|Hanl ll998: 



] l200l . 
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whether the formation of the second pulsar is consistent with a 
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